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For a 4-dimensional spatiaUy-flat Friedmann- Robertson- Walker universe with a scalar field (j){x), 
potential V{(j>) and constant equation of state w = p/p, we show that an expanding solution char- 
acterized by e = 3(1 -I- w)/2 produces the same scalar perturbations as a contracting solution with 
'e = 1/e. The same symmetry applies to both the dominant and subdominant scalar perturbation 
modes. This result admits a simple physical interpretation and generalizes to d spacetime dimensions 
if we define e = [{2d - 5) + (d - l)w]/{d - 2). 



I. INTRODUCTION 

In inflationary cosmology a nearly scale- invariant 
spectrum of primordial density perturbations is produced 
as comoving scales leave the Hubble horizon during an 
early burst of accelerated expansion 0] ■ In cyclic cosmol- 
ogy 3], the same perturbation spectrum is produced as 
comoving scales leave the Hubble horizon during a period 
of slow decelerated contraction 0, This agreement 
between two physically dissimilar models is unexpected, 
but not coincidental. As we shall show, the relationship 
between inflation and the cyclic model may be viewed as 
a special case of a surprisingly simple and general duality 
between expanding and contracting cosmologies. 

A Friedmann-Robertson- Walker (FRW) universe with 
a single scalar field 4> and potential V{(j)) is a simple 
yet important system. In particular, it is the canoni- 
cal 4d effective theory used to model the production of 
density perturbations in both inflationary and cyclic cos- 
mology. Recent results hint at a connection between two 
apparently-unrelated regimes of this model: (i) expand- 
ing, in the w —1 limit, and (ii) contracting, in the 
w —^ oo limit. Here w = p/ p denotes the ratio of pres- 
sure to energy density. Long-wavelength scale-invariant 
density perturbations are produced in the limits (i) and 
{ii) Q, and the small deviations from scale invariance 
near these two limits are related by the simple substitu- 
tion e ^ 1/e, where e = 3(1 w)/2 0. 

The relationship noted in 0] between expanding 
w ~ —1 models and contracting w ^ 1 models, turns out 
to be a special case of a general and exact duality relat- 
ing expanding and contracting models with identical per- 
turbation spectra. In this paper, we derive this duality, 
focusing on the case where e (or w) is time- independent. 
In a companion paper '^l, we generalize the discussion 
to the case where e is time- varying. When e is constant, 
or varies sufficiently slowly, the duality is simple: an ex- 
panding universe characterized by e produces exactly the 
same scalar perturbations as a contracting universe char- 
acterized hy t — 1/e. This duality applies in arbitrary 
spacetime dimension (not just 3-t-l); it applies for all w 



(not just the w —1 and w oo limits discussed in 
0, 0); it applies to all wavelengths (not just the long- 
wavelength limit); it applies to both the dominant scalar 
perturbation mode and a subdominant remainder (which 
are related to the growing and decaying modes, respec- 
tively) . 

This duality is of general theoretical interest since it 
provides a new relationship between expanding and con- 
tracting universes, and exposes an unexpected symmetry 
of cosmological perturbation theory. It is also relevant 
to cosmological models, like the ekpyrotic and cyclic Q 
scenarios, in which perturbations produced during a pe- 
riod of contraction are proposed to propagate through a 
bounce into a subsequent expanding phase. These mod- 
els require that the growing-mode long-wavelength per- 
turbation spectrum is preserved across the bounce. This 
has been a controversial matter. At first, some authors 
argued that growing-mode perturbations produced in a 
contracting phase must match to pure decaying-mode 
perturbations as one follows them across a bounce into 
an expanding phase 9]. At heart, this conclusion fol- 
lowed from requiring that the bounce corresponds to a 
comoving or constant-energy-density slice. However, re- 
cent five dimensional calculations 0, 0| indicate that 
comoving or constant-energy-density slices are inappro- 
priate for matching, since the bounce event (represented 
in five dimensions by a brane collision) is not synchronous 
in these slices. Matching on collision-synchronous slices 
results in the propagation of growing-mode perturbations 
across a bounce, with no change in the shape of the long- 
wavelength spectrum. Other aspects of the cyclic and 
ekpyrotic models have been criticized (see for 
replies), and some have argued that a bounce is impos- 
sible altogether While the consistency of a bounce 
remains to be proven, recent work has shown that the tra- 
ditional hazard of chaotic mixmaster behavior is strongly 
suppressed in the contracting phase of the cyclic model 
|14|. The metric perturbations exhibit ultralocal behav- 
ior in which anisotropics remain small, right up to a 
few Planck times before the bounce. In this situation, 
causality suggests that the bounce should not disturb 
correlations on macroscopic scales over which there can 
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be no communication in this finite time interval. Under 
these conditions, works by several groups sug- 
gest that perturbations generated during the contracting 
phase may pass into the expanding phase. The subject 
continues to be an area of active research. If the lat- 
ter suggestions are made rigorous, it would give added 
significance to the results presented herein. 

Other dualities have been identified in the literature 
[T^ HtL fl^ IT^ I20I I that relate cosmological solutions and 
perturbations. See Section IVll for a comparison. The 
duality presented here has the distinctive property that it 
relates two solutions that are stable under perturbations. 
Hence, both solutions can plausibly play a role in realistic 
cosmological models. 

The layout of the paper is as follows. In section 
we introduce the background (unperturbed) model: a 
spatially-flat FRW model with scalar field (j>, potential 
V{(j)), and constant e. In section Hill we briefly review 
scalar perturbation theory, before deriving exact solu- 
tions for Mukhanov's u and v variables. We note that 
u is invariant under e — > 1/e. This invariance is inde- 
pendent of our vacuum choice for the fluctuations, and 
provides our first glimpse of the duality. In section llVl we 
use u and v to separate scalar perturbations into pieces 
that are dominant and subdominant at long wavelengths, 
and show that each piece is both independently invariant 
under e —^ 1/e. We show how the dominant and sub- 
dominant pieces relate to the scalar perturbation grow- 
ing and decaying modes. In section^we consider tensor 
perturbations. In section IVII we contrast our duality 
with other kinds of cosmological dualities that have been 
studied in the literature El [Hill El 113 . In section 
IVIII we interpret our duality as a relation between the 
scale factor and the Hubble parameter, discuss its obser- 
vational significance, and mention some open questions. 
In an appendix, we generalize our results to d spacetime 
dimensions. 



II. BACKGROUND MODEL 

A spatially-flat Friedmann-Robertson- Walker (FRW) 
universe with scalar field cj) and potential Vicj)) is de- 
scribed by the metric 



ds2 = a(r)2 [-dr^ + dx^] 



(1) 



The unperturbed scalar field 4>a{T) and scale factor a(T) 
obey the Friedmann equations 



6(a'Va4) 
6(a'7a3) 



2p 

p-3p 



(2a) 
(2b) 



where we have chosen units such that c — h ^ 8ttG — 1, 
a prime (') denotes a conformal time derivative d/dr, 
and the energy density and pressure are given by 



p = (l/2)a-20^2 ^ y(<^o) 
p = (l/2)a-2^;,2 _ y^^^y 



(3a) 
(3b) 



Instead of the usual variable w 
convenient to use 

e = 3(1-1- w)/2. 



p/p, it will be more 



(4) 



to parameterize the equation of state. Equations H2I3() 
imply — 1 < w < 00 (or equivalently < e < 00). If w is 
near —1, then e <C 1 is the usual slow-roll parameter; but 
we make no slow-roll approximation in this paper, and e 
may be arbitrarily large. 

From now on, we shall assume that e is constant, not 
equal to unity. (The self-dual case e = 1 possesses spe- 
cial behavior which we shall not study here.) Then the 
solution of equations (|2I3|I is: 



a{r) 
Mr) 
U(0) 



±^^^ln|r| 



e - 1 
3 - e 

T^exp 



T(2e) 



1/2, 



(5a) 
(5b) 
(5c) 



where, to fix integration constants we have, without loss 
of generality, chosen the origin of conformal time so that 
a(0) — 0, normalized the scale factor so that a(l) = 1, 
and redefined <j)o ^ (/lo+constant so that 4>o{l) = 0. 
This solution separates into 4 cases: 

• (a) expanding, < e < 1, —00 < r < 0; 

• (b) expanding, 1 < e < 00, < r < 00; 

• (c) contracting, 0<e<l, 0<r<cx3; 

• (d) contracting, 1 < e < 00, —00 < t < 0. 

The Penrose diagrams for these 4 possibilities are shown 
in Figure^ Case (b) corresponds to an ordinary expand- 
ing FRW model with matter or radiation domination, 
e — 3/2 or 2, respectively. We are interested in the two 
cases, (a) and (d), in which t runs from —00 — ^ since, 
in these two cases, comoving length scales start inside 
the Hubble horizon at early times, and end outside the 
horizon at late times. Since we wish to study the am- 
plification of perturbations as modes leave the horizon, 
these are the two relevant cases. Thus, we will always 
assume e < 1 <^ expanding and e > 1 <^ contracting. 
The duality discussed below pairs solutions of type (a) 
with solutions of type (d). Figure ^ emphasizes that 
(a) and (d) have similar causal structure, but different 
singularities. 



III. SCALAR PERTURBATIONS 

In this section, we introduce relevant aspects of gauge- 
invariant scalar perturbation theory in 4 dimensions, and 
catch our first glimpse of the duality discussed in section 
IIVI For a thorough introduction to gauge-invariant per- 
turbations 01, see (23 and [l^. We work in Fourier 
space throughout, so every perturbation variable carries 
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where fc = | fc [ is the magnitude of the (comoving) Fourier 
3- vector. On the other hand, the gauge-invariant pertur- 
bation variable ( is most easily understood in "comoving 
gauge" {Ht = 5T^i = 0), where it represents the cur- 
vature perturbation on spatial-hypersurfaces, and is re- 
lated to the spatial metric perturbation in a simple way: 
C = —Hl. The condition 5T^^ = also implies that 
(5(/) = in this gauge. C, obeys the equation of motion 

C" + 2{z'/z)C' + k'C^O (8) 

where z = a^^p/a'. $ and C are related to each other by 

C = $+i[(a/a')^' + *] (9a) 
$ = -e{a' /a)k~^C' ■ (9b) 



FIG. 1; Penrose diagrams for spatially-fiat FRW universes 
with: (a) < e < 1, expanding; (b) 1 < e < oo, expanding; 
(c) < e < 1, contracting; (d) 1 < e < oo, contracting. 
The left edge of each diagram is the world line of a comoving 
observer at the origin; curved lines represent other comoving 
world lines and spatial hypersurfaces. The Hubble horizon is 
a curve connecting the 90° vertex to the lightlike boundary, 
but the precise curve depends on e. For illustration, we have 
shown the horizon for e = in (a, c) and for e = 2 in (b, d). In 
this paper, we focus on cases (a) and (d), in which comoving 
scales exit the Hubble horizon. 



an implicit subscript k which, for brevity, is not shown 
explicitly. Write the perturbed metric 



(1 + 2AY)At^ - 2BY,dTdx' 
- [(1 + 2HlY)5,j + 2HtY,j] dx'dx^ (6a) 



and perturbed scalar field 

<P = 0o(t) + > 



{r)Y 



(6b) 



where Y{x), l^faf), and Yij{x) are scalar harmonics (see 
Appendix C in |23|'). The corresponding perturbations 
of the Einstein tensor and energy-momentum tensor (see 
Appendices D and F in are related to one another 
through the perturbed Einstein equations, 5G^^ — ST'^^. 

It is well known that scalar perturbations in a spatially- 
flat FRW universe with scalar field (j) and potential V{(f)) 
are completely characterized by a single gauge-invariant 
variable. But the choice of this variable is neither unique 
nor standard; two of the most familiar options are the 
"Newtonian potential," and the "curvature perturba- 
tion," C- 

The gauge-invariant Newtonian potential $ is most 
easily understood in "Newtonian gauge" {B = Ht = 0), 
where it is related to the metric perturbations in a simple 
way: ^ — A = —H^. It obeys the equation of motion 



$" + 2 




$' + 2 




a 00 





2H' - 2n^ 



$ = (7) 



Note that our definitions for $ and C agree with those 
in ,23]. But beware: in 22], the gauge-invariant Newto- 
nian potential is denoted 5*, while $ denotes a different 
(though closely related) variable. 

It is convenient to introduce new variables, u and v 
|23 ]. by multiplying <I> and C by fc- independent functions 
of r: 



u = (a/0o)$ V = z(. 



(10) 



Note that u and v have the same fc-dependence as $ and 
C, respectively, and may serve as surrogates for $ and ^. 
If we also define background quantities, 9 and z: 



1/Z: 



then u and v obey simple equations of motion 



u 
v" 



{k^ - e"/e)u 

(fc2 - z"/z)v 



and are related to each other by 



kv = 2k[u' + iz'/z)u] 
1 

2k' 



-ku 



{e'/6)v] 



(11) 



(12a) 
(12b) 



(13a) 
(13b) 



We must choose a vacuum state for the fluctuations, 
which corresponds to specifying appropriate boundary 
conditions for u and v (see Ch.3 in Birrel&Davies [2^). 
The standard choice is the Minkowski vacuum of a co- 
moving observer in the far past (when all comoving scales 
were far inside the Hubble horizon) , corresponding to the 
boundary conditions 



u 

V 



i(2fc)-3/2e- 
{2k)-^/\- 



■ikr 



(14a) 
(14b) 



as T — oo. Using (|13|l . it is easy to check that these 
two boundary conditions are equivalent. 
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When e is time-independent, we can use (0) to find 



z"lz 



(e-l)2r2 
2 - e 



15a) 
15b) 



(e- l)2r2 
Then we may solve (I12f) to obtain 

M(a;) = (x) + ^(2)^(2) (^^1 (^g^^ 

where x = fc|T| is a dimensionless time variable, yl''^'2) 
and i3(^'2) ^re constants, H^'^\x) are Hankel functions, 
and we have defined 



a= ^/{9"/9)t^ + 1/4 



/3= 



e + 1 



e - 1 
e-3 



1 



(17a) 
(17b) 



In the far past (x — > 00) we use the asymptotic Hankel 
expression. 



so the boundary conditions (|14|l imply 

u = g(^x/4fc)i/2ff(i)(2;) 
v = V2{^x/Ak)^'^H^^\x) 



where 



Vi = exp[i(2a + 3)7r/4] 
=exp[i(2/3 + l)7r/4] 



(18) 

(19a) 
(19b) 



(20a) 
(20b) 



are fc-independent complex phase factors. 

Note from (|15a|l that the equation of motion for m, 
(|12a|l . is invariant under e — > 1/e, while the boundary 
condition, H14a(l . is independent of e. As a result, our 
expressions ()17a|l for a and (|19a|l for u are invariant under 
e — > 1/e. This is our first glimpse of the duality discussed 
below. 

We stress that this result does not depend on the par- 
ticular vacuum choice (|14a|l . Any boundary condition 
that is independent of e (or, more generally, invariant 
under e — > 1/e) will work. And it is natural to expect 
the boundary condition to be independent of e, since it 
is imposed in the far past, when comoving scales are far 
inside the Hubble horizon. 



IV. DOMINANT AND SUBDOMINANT MODES 

In this section, we show that u and v can be decom- 
posed into pieces that are dominant and subdominant 



at long wavelengths such that each is invariant under the 
transformation e — *■ 1/e. The dominant and subdominant 
parts are closely related growing and decaying modes over 
the range of w relevant to cosmological model-building, 
as we explain below. 

For this analysis, it is convenient to scale u and v by 
appropriate powers of |r|, so that they only depend on k 
and r through the dimensionless combination x = k\T\. 
Thus, using (|19al I19b|l . define 



(21a) 



u^W'u = g(V4)V2ij(i)(,) 

v=\t\-^'''v = V2{^/^Y'^Hf\x). (21b) 

Note that u and v have the same fc-dependence as u and 
w, respectively, or and respectively, and may be used 
in place of these more standard variables. 

To make the meaning of "dominant" and "subdom- 
inant" precise, consider two linearly independent func- 
tions fi{x) and f2{x). li\m\x^o f2{x) / fi{x) exists, then 
/i (x) and /2 (x) can be related by a linear transformation 
to two new functions, /dom(a;) and fsuhix), satisfying 



lim /sub(a;)//dom(a;) 

a:— ^0 



0. 



(22) 



So the subdominant piece, /sub(a;), becomes negligible 
relative to the dominant piece, /dom(a;), for small x {i.e. 
far outside the horizon). The condition (|22|l uniquely 
determines the subdominant piece (up to an overall nor- 
malization constant) to be 



/sub (a;) = f2{x) - 



h{y) 



hm 

a^" fi{v) 



(23) 



but does not uniquely fix the dominant piece. Rather, 
/dom(a;) may be any linear combination of fi{x) and 
f2{x) that is hnearly independent of fsuh{x). We can 
now choose fi{x) — u{x) and f2{x) — v{x), and find the 
corresponding dominant and subdominant functions. 

Let us first calculate fsuh{x). To compute 
linia^^o /2 (a;) //i (x) , recall the Hankel identity 



HP (a 



-*[r(s)/7r](x/2)- 



0, (24) 



where s > and r(s) is the Euler gamma function. Then 
from (|21al21b|l we find 



lim 



/2(X) 



hix) 



{V2/Vi)Aa 




< e < 1 

1 < e < 00 



(25) 



where we have used the fact that a + 1 = (3 when e < 1 , 
while a + \ > (3 when e > 1. Now, substituting (|21a|l . 
(|2rbp and into (^5)) , using the Hankel identities 



hP,{x)+h'-1^,{x)^^hP{x) 



h[^}{x) =e''^^ijy^(x) 



(26a) 
(26b) 



and paying careful attention to absolute value signs, we 
find 



/sub(a;)=7'3(7r/4)i/2i/W(x) 



(27) 
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where we have defined 



7 



la -11 



(28) 



and Vj, = exp[i(27 + l)7r/4] is a fc-independent complex 
phase factor. 

Now let us turn to /dom(a;)- The fact that 
linia^^o fsnh(x) lv(x) — when e < 1 shows that the dom- 
inant mode contributes to v in an expanding universe. 
Since /sub (2;) = v{x) when e > 1, w is purely subdomi- 
nant (i.e. contains no dominant-mode contribution) in a 
contracting universe. By contrast, u and /s(x) are always 
linearly independent, and hence \mix—,Q fsuh{x)/u{x) = 
for all e. Thus, we can use the freedom in defining /dom 
to choose 

/dom(x) = U{X) = g(^/4)l/2ij(l)(:,) (29) 

for the dominant piece. 

Notice that the expressions (|27|) for fsuh{x) and H29|) 
for /dom (2;) are both invariant under e — > 1/e, because a 
is invariant. Thus, we see that the subdominant mode 
automatically possesses this symmetry, since /sub (2;) is 
uniquely determined (up to a normalization factor) by 
the condition H22() . Furthermore, we have shown that we 
can choose a linear combination of u{x) and v{x) such 
that /dom (a;) displays the same symmetry. (If we had 
made the wrong choice for /dom (2;), then the exact e 
1/e symmetry of the dominant mode would be hidden, 
and would only re-appear in the long- wavelength limit.) 

In the a; — > limit, define dominant and subdominant 
scalar spectral indices, ridom and risub, which satisfy 



.ml' OC 

a^^l/subP oc x" 



a^^l/domT X" 



(30a) 
(30b) 



Using |(57l and (|21l , along with the Hankel identity (|23| , 
we find 



l = l-2a = 1 



1 



"sub -1=3 - 27= 3- 



e - 1 
e+ 1 



1 



- 2 



(31a) 
(31b) 



These spectral indices are plotted in Fig. (as a func- 
tion of w), and in Fig. 13) (as a function of e). Again, no- 
tice that ridom and risub are both invariant under e 1/e. 
This symmetry is manifest in Fig. ^p. 

Since e lies in the range < e < 00, this duality for- 
mally pairs every expanding (e < 1) universe with a con- 
tracting (e~ > 1) universe, and vice versa. However, the 
background solution (|3J is only stable against small per- 
turbations in two cases: (i) expanding with e < 1 or 
{ii) contracting with e > 3 ,6i, 14J. Thus, an expand- 
ing model and its contracting dual are both stable when 
e < 1/3 {w < -7/9) and ? > 3 (w > 1). Also note, 
in agreement with |^, that there are only two limits in 



which an approximately scale-invariant (rig 



1 



0) 



□3-1 



subdominant mode 
(expanding) 



(a) 




subdominant mode 
(contracting) 



W 



dominant mode 
(contracting) 



dominant mode 
(expanding) 



subdominant mode 
(contracting) 



n - 1 



subdominant mode 
^^1^ (expanding) 



(b) 



dominant mode 
(contracting) 




dominant mode 
(expanding) 



FIG. 2: The dominant and subdominant scalar spectral in- 
dices (a) as a function of w and (b) as a function of Ine. Note 
especially the symmetry of (b). 



e — > {w — > —1), corresponding to the inflationary 
regime and (ii) when e — s- 00 (w — s- 00) corresponding 
to the cyclic/ekpyrotic regime. 

The dominant and subdominant pieces, /dom (2^) and 
/sub (a;) are related to the growing and decaying modes 
of u and v, which may simply be obtained by replacing 
the Hankel functions in I|21al21b(l with the corresponding 
Neumann functions or Bessel hmctions: 



^grow (a-) 
Wgrow OC Ys(x) 



^dccay ^ ^ <-^a(*^) 
^'dccay (X Jf3{x). 



(32a) 
(32b) 



spectrum of scalar perturbations is produced: (i) when 



We may define growing-mode and decaying-mode spec- 
tral indices for u and v, just as we did for /sub and /dom 
in Now restrict attention to the ranges e < 1/3 and 
t> 3 {i.e. the range over which the duality pairs stable 
expanding models to stable contracting models). Then, 
using H32I) along with the (|24|l . we find 

• The growing-mode and decaying-mode spectral in- 
dices associated with u are simply equal to ndom 
and risub, respectively. 

• The growing-mode spectral index associated with 
V is equal to ndom in an expanding (e < 1/3) uni- 
verse, but is equal to risub in a contracting (e~ > 3) 
universe. 
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V. TENSOR PERTURBATIONS 

Tensor perturbations are much simpler than scalar per- 
turbations. The perturbed metric is 

dsVa^ = -dr^ + + 2hTY}^^]dx'dx^ (33) 

where Y^^^ is a tensor harmonic (see Appendix C in |22||). 
The tensor perturbation hr is gauge-invariant and obeys 



hT + 2{a'/a)h'j, + k'^hT = 0. (34) 
It is useful to define a new variable /t = ahx which obeys 



+ (P - a"/a)fT = 0. 



(35) 



Again, the standard vacuum choice is the Minkowski vac- 
uum of a comoving observer in the far past, correspond- 
ing to the boundary condition 



/t - (2fc) 



-l/2g-ifcT 



-OO. 



(36) 



We can now solve for fr, just as in the scalar case. But it 
is simpler to notice that lO and imply z(t) oc a(r), 
and hence z" / z — a" /a when e is constant. Thus, 
since v and /t obey identical equations of motion (com- 
pare (|12b|) with (123) and boundary conditions (compare 
((T4E|l with (ISnj), we find 



(37) 



The tensor spectral index is defined in the long- 
wavelength limit by fc^l/Tp oc fc"^. Using and ||1T7|) 
we find 



riT = 3 - 2/3 = 3 



e - 1 



(38) 



Note, in particular, that this expression is not invari- 
ant under e ^ I/e. An expanding universe with equation 
of state e produces a tensor spectrum which is much red- 
der than the tensor spectrum produced in a contracting 
universe with e~= 1/e: ut < nx — 2. 



VI. OTHER DUALITIES 

It is interesting to contrast our duality with other cos- 
mological dualities that have been discussed in the liter- 
ature. 

One duality, due to Wands 0| (see also P^irs 
models that share the same "w" perturbations. By con- 
trast, our duality pairs models that share the same "u" 
perturbations. (Note: the variable called "it" in is 
called "z;" in o ur p aper, in agreement with Mukhanov's 
convention |23,I3|-) Thus, whereas our duality connects 
expanding and contracting models through the substi- 
tution e — > ? = 1/e (which leaves 0"/6, and hence u, 



invariant). Wands' duality instead uses the substitution 
e — > ? = (2e — 3)/(e — 2) (which leaves z"/z, and hence v, 
invariant). For example, his duality pairs an expanding 
inflationary solution (e = 0, w ~ —1) with a contract- 
ing dustlike solution (e ~ 3/2, w — 0). Another differ- 
ence between our duality and Wands' stems from the fact 
that V is purely subdominant (i.e. contains no dominant- 
mode contribution) in a contracting universe (see section 
IIV|) . Thus, whereas our duality maps the expanding- 
phase dominant mode to the contracting-phase dominant 
mode, and the expanding-phase subdominant mode to 
the contracting-phase subdominant mode. Wands' du- 
ality instead associates the expanding-phase dominant 
mode with the contracting-phase subdominant mode. 

A second interesting duality, discussed by Brustein et 
al. applies to a broad class of cosmological pertur- 
bations. Associated with each type of perturbation is a 
"pump" — a particular function of the background fields. 
The Hamiltonian governing a given perturbation is in- 
variant under a transformation that swaps the perturba- 
tion with its conjugate momentum, and simultaneously 
inverts the appropriate pump function JlS] ■ It is instruc- 
tive to apply this duality to the u and v variables consid- 
ered in this paper. In this case, the Brustein et al. duality 
associates an expanding solution characterized by 6 (or 
e) to a contracting universe characterized by 6* = 1/9 — z 
(or ? = 2 — e) . The transformation e ^ t — 2 — e effec- 
tively swaps the variables u and v 



u —> u = {i/2k)v 
V ^ V ^ (2k/i)u 



(39a) 
(39b) 



as may be verified from (|19all9b|l . 

Recall that e lies in the range < e < oo. Thus, 
our duality formally pairs every expanding (0 < e < 1) 
solution with a contracting (1 < e~ < oo) solution, and 
vice versa. By contrast. Wands' duality relates every 
expanding solution to a contracting solution with 1 < 
t< 3/2; but contracting solutions with 3/2 have no 
expanding dual. Similarly, Brustein et al.'s duality pairs 
every expanding solution with a contracting solution in 
the range 1 < e~< 2; but contracting solutions with 2 
have no expanding dual. 

The constant-e background solutions lO are only prac- 
tically useful if they are dynamically stable. Recall 
that the contracting solutions are stable if and only if 
e > 3 (w > 1) la Ell. Thus, Wands' and Brustein 
et al.'s dualities relate every expanding solution to an 
unstable contracting solution. By contrast, our dual- 
ity relates stable expanding solutions with < e < 1/3 
(—1 < w < —7/9) to stable contracting solutions with 
3<e~<oo(l<w< oo). In terms of the spectral 
index, this means that Us > may be produced either 
by a stable expanding phase or by a stable contracting 
phase. In the real universe, the condition tIs > is sat- 
isfied (experiments favor Ug « 1), so that our duality 
is of practical relevance in cosmological model building. 
Some properties of the three different duality relations 
are summarized in Table I. 
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Transformation of 


Range 


Maps stable 




background 


perturbations 


or e 


to stable? 


Ket. 




u u 


[0,00) 


Yes, if Jis > 




2e-3 

6-2 


i; — > i; 


[0,3/2] 


never 


m, 17] 


e ^ 2 - e 


u « (i/2k)v 


[0,2] 


never 


fl8] 



where t is the proper time of a comoving observer and 
iJ-Ms the "Hubble scale." Thus, two dual universes are 
related by 



TABLE I: Comparison of the duality presented here with 
those presented by Wands [iMllTl and by Brustein et al. 
The first two columns show how the background and per- 
turbation variables transform under each duality. The third 
column shows the range of e to which the duality applies. 
The fourth column indicates the condition under which an 
expanding background solution and its contracting dual are 
both stable under small perturbations. 



Finally, a number of authors have discussed cosmolog- 
ical symmetries of the low-energy string effective action. 
If one neglects all fields in this action besides the dilaton 
and the metric, then there is a well-known "scale-factor 
duality" starting with any cosmological solution, 

one can use this duality to generate new solutions. If, 
in addition to the dilaton, one includes other fields (ax- 
ions, moduli,...), then the cosmological solutions may 
display more general dualities, and the resulting pertur- 
bation spectra may be invariant under these dualities 
[2(1] . But note that these symmetries typically relate dif- 
ferent solutions of a single effective action. By contrast, 
the dualities in Table I relate two different cosmological 
background solutions corresponding to two different La- 
grangians: an expanding universe, with potential V{(f>), 
is dual to a contracting universe, with a different poten- 
tial V{cj)). 



VII. DISCUSSION 

Beyond its inherent theoretical interest, our duality 
may be observationally relevant if (as discussed in section 
m long-wavelength correlations produced during a con- 
tracting phase successfully propagate into a subsequent 
expanding phase. This suggests a fundamental degener- 
acy: an ideal measurement of the "primordial" scalar per- 
turbation spectrum may be unable to determine whether 
the perturbations were generated by an expanding phase 
or by its contracting dual. Luckily, as shown in sectionlVI 
tensor perturbations break this degeneracy: an expand- 
ing model produces a much redder tensor spectrum than 
its contracting dual. In particular, a detection of ten- 
sors in the cosmic microwave background would indicate 
that these perturbations were generated in an expanding 
phase, since the dual contracting phase would produce 
an undetectably small tensor spectrum on these cosmo- 
logical length scales |26| . 

Using the background solutions ||SJ), we may think of 
the duality as relating two different scale factors, a{T) and 
a(T), or two different scalar potentials, V{4)) and V{(1)). 
Alternatively, recall that a oc \t\^^^ and = a? /a' oc 



d XnH/d Ina = d Ina/d XnH . 



(40) 



So the duality effectively swaps the scale factor with 
the Hubble scale, and simultaneously swaps expansion 
and contraction. For example, in the e ^ (ui — > — 1) 
limit, the scale factor grows rapidly while the Hubble 
length grows slowly; whereas in the e oo (w — > oo) 
limit, the Hubble length shrinks rapidly while the scale 
factor shrinks slowly. Expanding models in which modes 
exit the horizon most rapidly {w — > —1) and most slowly 
(w —>■ —1/3), are associated with contracting models in 
which modes exit most rapidly (w — > oo) and most slowly 
(w —1/3), respectively. 

Finally, let us mention several open questions. We have 
studied the duality using a simple model — spatially-flat 
FRW spacetime with a single scalar field (p; what happens 
in more complicated models? We have used linear per- 
turbation theory; what happens in the nonlinear regime? 
We have restricted ourselves to time-independent e; what 
if we allow time- varying e? This final question is treated 
in a companion paper [g] . 
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APPENDIX: GENERALIZATION TO 
ARBITRARY SPACETIME DIMENSION 



In sections iHl— fVl we restricted our discussion to 4 
spacetime dimensions. Here we sketch the generalization 
to d spacetime dimensions, for d > 4. To understand 
this appendix, read sections Ull— TVl first . Gauge-invariant 
perturbation theory in d spacetime dimensions is treated 
thoroughly in 22] (see especially the appendices therein). 

The background metric obeys the Friedmann equa- 
tions 

{d-l)id-2)a'ya'^ = 2p (A.la) 
{d -l)id- 2) a"/a^ = (5 - d)p - {d - l)p (A.lb) 

where p and p are given by Instead oi w = p/p, 
parameterize the equation of state with 



(2(i- 5) + (d- l)ui 
d-2 ■ 



(A.2) 
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From eqs. (|51 lA.l|l we find —1<w<ooot j-^ < e < oo. 
For constant e, the solution of lA.ip is 



l(r) = |r|2/[(rf-2)(e-l)] 
1 




(A.3a) 
(A.3b) 

(A.3c) 



wfiere we fiave cliosen a(0) = 0, a(l) = 1 and 4>o{l) = 0. 

As sliown in |22j , perturbations in d spacetime dimen- 
sions may be decomposed into scalars, vectors and ten- 
sors, just as in 4 dimensions, and gauge-invariant vari- 
ables may be defined. In particular, we can again in- 
troduce scalar perturbations through equations I|t)all6bp . 
and describe these perturbations with a single gauge- 
invariant variable. The gauge-invariant Newtonian po- 
tential, <&, is most easily understood in "Newtonian 
gauge" {B = H]^ = 0), where it is related to the met- 
ric perturbations in a simple way: ^ ~ A~ — (d — 3)i?i. 
The gauge-invariant variable C is most easily understood 
in "comoving gauge" {Ht = 5Tf = 0), where it is re- 
lated to the spatial metric perturbation in a simple way: 
C, = — i?L- Note also that 50 = in this gauge. If we 
take the d-dimensional definitions of u, v, 9 and z to be: 



I fd-2\ a('^-2)/2 

u = — I $ 

2\d-3j cj)'. 



l/z 



v = zC (A.4) 



(A.5) 



then sections IIIII and ll VI immediatelv generalize d dimen- 
sions. Indeed, it is straightforward to show that equa- 
tions jnil through all remain true in d dimensions, 
provided we take e, u, v, 9, and z to be given by their 
d-dimensional definitions (jA.2|l . (jA.4p and (jA.5|l . 



In particular, this means that the duality extends to 
d-dimensions: u, /grow and /decay are all invariant under 
e 1/e. The duality pairs expanding solutions in the 
range < e < 1 with contracting solutions in the 
range 1 < e < cxd. Note that for d > 4, there are some 
contracting models with no expanding dual; but for d = 
4, every expanding model has a contracting dual, and 
vice versa. 

In the a; — > limit in d dimensions, define the growing 
and decaying spectral indices 



I /grow I CC X ^ 



f grow 

I /decay | CK X ^ 

Then, using lfT7a|) . lEH), 117|), 1(13 and (EHl) we find 



(A.6a) 
(A.6b) 



(A.7a) 

"decay - i = d-1- ^—^ - 2 (A.7b) 
■ e — 1 

with e given by (|A.2p . Note that these expressions are 
invariant under e — » 1/e. 

As shown in p3 |. we may again introduce tensor per- 
turbations through eq. (|33|l . Now define the gauge- 



invariant variable /^ = a''* ^^/^/it which obeys 



(A.8) 



along with the boundary condition H36() . Using H36|l . 
(|A.3p and (|A.8p we find that solution for /t is still given 
by H37|l in d dimensions. In the k limit, define the 
tensor spectral index /c'^^^l/P oc fc"^. Then, using H24|) 
and lEZl we find 



riT 



1 



(A.9) 



with e given by l)A.2(l . Note that this expression is not in- 
variant under e — > 1/e. This completes the generalization 
to d spacetime dimensions. 
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